Optical properties of an effective one-band Hubbard model for the cuprates 
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We study the Cu and O spectral density of states and the optical conductivity of Cu02 planes using 
an effective generalized one-band Hubbard model derived from the extended three-band Hubbard 
model. We solve exactly a square cluster of 10 unit cells and average the results over all possible 
boundary conditions, what leads to smooth functions of frequency. Upon doping, the Fermi energy 
jumps to Zhang-Rice states which are connected to the rest of the valence band ( in contrast to 
an isolated new band in the middle of the gap). The transfer of spectral weight depends on the 
parameters of the original three-band model not only through the one-band effective parameters but 
also through the relevant matrix elements. We discuss the evolution of the gap upon doping. The 
optical conductivity of the doped system shows a mid-infrared peak due to intraband transitions, 
a pseudogap and a high frequency part related to interband transitions. Its shape and integrated 
weight up to a given frequency ( including the Drude weight ) agree qualitatively with experiments 
in the cuprates for low to moderate doping levels, but significant deviations exist for doping x > 0.3. 



PACS numbers: 74.72.-h, 78.20.-e, 79.60.-i, 75.40.Mg 



I. INTRODUCTION 

In recent years, there has been much interest in the 
spectral and related electronic properties of cuprates su- 
perconductors and similar materials p]- |31j| . In spite of 
the considerable research effort, there are several issues 
which remain to be clarified. For example, photoemission 
experiments at optimal doping show that electrons 
have a large Fermi surface of area (1 — x) where x is the 
amount of hole doping, while at small doping it is ex- 
pected that the Fermi surface consists of four small hole 
pockets centered at (±7r/2, ±7r/2) of total area x. The 
evolution of the Fermi surface remains a tough problem 
[ 31 1 . There is a recent theoretical study on this subject 

Another related issue is the appearance and evolu- 
tion for i ^ of states at energies which lie in the 
gap for x = 0. There are at least two physical pic- 
tures. On the basis of several measurements, and par- 
ticularly the change of sign of the Hall constant with x 
in I^-ajSr^CuO^ Sreedhar and Ganguly || proposed 
that the gap in the spectral density gradually closes with 
increasing x, and after the closing, the system has com- 
pleted its evolution from a Mott ( or better said "charge- 
transfer" ) insulator to a band metal. As soon as doping 
begins the Fermi level jumps to the valence band. This 
jump is consistent with the measured x-ray absorption 
spectrum (XAS) |||7]]. Instead, from the optical conduc- 
tivity H and high-energy spectroscopies ||, it has been 
suggested that a new band of midgap states is formed, 
taking spectral weight from both, conduction and va- 
lence bands, and the Fermi level remains in the midgap 
band. A minimal shift of the Fermi level with doping 
has been observed by photoemission in L^-zSr^CuO^ 
H and in the electron doped Nd2- :E Ce a: Cu04 M, but 



not in F^S^Cai-zYzCuO^Cs |l0[ . It might be possi- 
ble that the states in the middle of the gap are created by 
deep donor levels originated by the substitution of La +3 
by Sr +2 ions. The resulting "impurity-state model" is 
consistent with the experimental evidence in the oxide 
(Nd,Sr)Co0 3 Q. The effect of this substitution has 
been discussed in Ref. 0. However, inhomogeneities are 
usually not included in the theoretical treatments, and it 
seems necessary to generalize the Hubbard model to ex- 
plain the electronic structure of some non-stoichiometric 
oxides @0. 

To discuss the validity of the translationally invari- 



ant three-band Hubbard model 1 33 1 or effective models 
derived from it p9| , p0| , p4| -|3r|| in the description of the 
cuprates, it is necessary to know precisely the spectral 
properties of these models. One of the most used and ef- 
fective ways to study these properties in two dimensions 
at zero t emp erature is the exact diagonalization of small 
clusters |i^pj|-^2]j24|j29|l . Among the different effective 
models, those similar to the t — J one 
have the smallest number of states per unit cell. How- 
ever, like the spin-fermion ( or Kondo-Heisenberg) model 
|]l9| , po| , their Hilbert space is too small to allow a de- 
scription of both, valence and conduction bands. In con- 
trast, the number of states per unit cell of the three-band 
model is so large that the largest exactly solved cluster 
contains only four unit-cells ]l7],[l8|,|3| . Thus, the most 
apropriate effective model to study the evolution of the 
gap seems to be the extended Hubbard one J35|. Nu- 
merical studies of the Hubbard model show a change of 
sign of the Hall constant as a function of doping, but it 
has not been related with a reconstruction of the spectral 
density |37|-p0|] . The spectral properties of the Hubbard 
model have been calculated in periodic square clusters of 
8,10 and 16 sites |2U-E3]. Without artificial broadening, 
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the result of these calculations is given by a set of delta 
functions at different frequencies ( see for example Fig. 11 
of Ref ^2|), from which it is not possible to distinguish 
if in the thermodynamic limit there is one, two or more 
gaps in the spectral density of states. Also the Drude 
weight shows important finite-size effects, being negative 
near half-filling. 

In this work we solve exactly the effective one-band 
extended Hubbard hamiltonian in a square cluster of 10 
unit cells and average over all possible boundary condi- 
tions Q| . This technique has been also used by Poilblanc 
to study the optical conductivity of the t — J model |24| ] 
and by Xiang and Wheatly |^2| to calculate the disper- 
sion relation of one hole in a generalized t — J model, ob- 
taining good agreement with results of the self-consistent 
Born approximation and experiments in Si^CuC^C^. To 
our knowledge this approach has not been applied to 
the spectral properties of Hubbard or extended Hubbard 
clusters. It allows us to obtain continuous spectral densi- 
ties and optical conductivities and according to our stud- 
ies in small rings ( see section III ), the convergence to 
the thermodynamic limit is faster. 

It must also be pointed out that the spectral proper- 
ties of the Hubbard model derived from the numerical 
methods mentioned above, as well as analytical studies 
of the change in the spectral weights with doping ||(]] 
do not directly correspond to the cuprates as described 
by the original three-band Hubbard model. To calculate 
any property of this model using an effective Hamilto- 
nian, one should transform the corresponding operators 
from the three-band to the effective Hamiltonian. This 
has been done for the spin- fermion @, Hubbard @ and 
t — t' — t"—J |3(| as effective models for the cuprates. In 
particular, Feiner explained the electron-hole asymmetry 
in the spectral weight in the cuprates ( which is absent in 
the " isolated" Hubbard model) jij] . Here we use and dis- 
cuss the transformed operators for spectral density and 
optical conductivity. 

In section II, we briefly review the extended one-band 
Hubbard model as an effective model for cuprate super- 
conductors, and construct the relevant operators. In sec- 
tion III, we discuss the numerical techniques and the 
method of integrating over arbitrary boundary condi- 
tions. In section IV, we present the results for the spec- 
tral density, while section V contains the results for the 
optical conductivity, Drude weight and related quanti- 
ties. The conclusions are presented in Section VI. 



II. THE EFFECTIVE EXTENDED ONE-BAND 
HUBBARD MODEL AND TRANSFORMED 
OPERATORS 

Experimental evidence about the symmetry of holes in 
high-T c superconductors |3 45 4C ] , as well as constrained- 
density- functional calculations [4^,Q, support the ap- 



propriateness of the three-band Hubbard model |33| for 
the description of these systems. The Hamiltonian is: 



H 3b = AJ2p]*Pj<? + tpd^2(pt +Sa div + h.c.) - (1) 

jo iSa 

iSocr' 



where d ia (Pj a ) creates a hole in the d x 2_ y 2 ( p a ) or- 
bital at site i ( j ) with spin c, i + S ( j + 7 ) label 
the four O atoms nearest-neighbors to Cu (O) site i ( 
j ). The phases of half the d and p orbitals have been 
modified by a factor (-1) in order to have t p d > and 
t pp > independently of the direction. While other Cu 
and O orbitals should be included to explain Raman ex- 
periments p9| , ^p| , these additional states affect neither 
low-energy spectral properties (below the bottom of the 
p n band |5C ] ) nor the validity of the one-band effective 
model Hu |fj. 

The transformation of the low-energy part of H 3b to 
the effective model Hu, can be summarized in four steps 
p4] , j35 50, 5|]: (i) change of basis of the O orbitals p*,p V j 
to orthogonal Wannier functions a>i,ji centered at each 
Cu site i, with symmetries b\ g ( the same as the d x i_ y i 
orbital ) and a\ g respectively; (ii) exact solution of the 
cell Hamiltonian Hi ( the terms of H$b which contain 
only operators acting on cell i ); (hi) mapping the low- 
energy states I to) of H^b into corresponding ones |to) of 
the Hilbert space of the one band model (|m) <-> |m) 
); (iv) transform the operators Hi and H 3b — to 
the new basis |to). A transformed operator O can be 
expressed in terms of the matrix elements of the original 
operator O as: 



O = J2(n\0\m)\n)(fh\ 



(2) 



Finally, the effect of states of H^ which do not have cor- 
responding ones in Hu, can be included as perturbative 
corrections |3q| , but we will not do it here. In Fourier 
space, the p a orbitals are related with the O Wannier 
functions by: 



y 



/3 k (cOS (-^— )a k a - COS (— 7k<r) 
= /3 k (cOS (-y)a kj + COS (-^|— ) 7ka) 



with /3 k 



(1 + i cos (k x a) + i cos (k y a)) 



(3) 



The 



eigenstates of Hi retained in the mapping procedure are ( 
in addition to the vacuum at site i ), the lowest one-hole 
doublet \ia) and the lowest two-hole doublet \i2). They 
can be expressed as: 



2 



i2) = [^=(< 



A^yim (4) 



where 



[Bidt-^aJjO) 



These states are mapped into those of the one-band 
model by the correspondence \ia) <-> c\ a \0), \i2) <-> 
cj^ctjO). This assignment and Eq.(2) applied to 
leads to the effective one-band generalized Hubbard 
Hamiltonian plfl , with occupation dependent nearest- 
neighbor hopping and repulsion ( interactions and hop- 
pings at larger distances are neglected for simplicity ): 



Di 
D 2 
Pi 
P 2 



(iO\d l<7 \ia) = B x 

(ia\d la \i2) =(^L+A 3 B 1 ) 

(iO\a ia \i<j) = -B 2 

,A X B X 



(8) 



(ia\a irT \i2) = -(- 



V2 



A 2 B 2 ) 



Hib = Ei n ia + U n^riii + ^ c f ja c ia {t AA (l - ni, a )(l - n^ a ) + t AB [ni, 9 (l - n^ a ) + n jy9 (l - n ha )} + (5) 



<ij>o 



tBBn ttS nj, s } + h.c. + {V11 (1 - n uS ) (1 - n^ s i) + V%%n i>a nj^i 

Vi 2 [(l ~ rii.^n^g, + n i:< j(l - n^/)]}^/^ 



Several particular cases of this model have been widely 



studied |2|-§8|. For t AA + t BB = 2t AB and Vy = 0, the 
model is shown to lead to superconductivity for small 
U and adequate doping [52]| , particularly in one dimen- 
sion l5^-|5^1- For Vij — t AB — 0, the model has been 
exactly solved in one dimension (56). For t AB = 0, 
V 22 = V21 = V11 and half-filling, the exact ground state 
has been found in a wide range of parameters for arbi- 
trary lattices in arbitrary dimension, and the phase di- 
agram separating metallic, Mott insulating and charge- 
density waves regions has been established [ |57| . 

The spectral properties in two dimensions have been 
studied numerically in the charge-density-wave ( large Vij 
) regime |Q , and for Vij = and small t AB at half- filling 
p9[ . The model has been also used recently to study the 
regions of phase separation and valence instabilities of 
the original three-band model . Typical values of the 
parameters of and Hu are given in tables I and II 
respectively. 

At zero temperature, the Cu and O spectral densities 
are defined by: 



po(w) = ^\(n\Q k \g)\*8(u- 

k.n 



E„ 



(6) 



where \n),E n are the eigenstates and energies of Hn,, \g) 
labels the ground state for a given number of particles TV 
and L is the number of unit cells in the system. For elec- 
tron spectral densities of Cu or one of both O atoms in 
the unit cell ( measured by inverse photoemission ) , Ok is 
the operator in the representation of Hu of dk or p x k ' v re- 
spectively. For hole spectral densities ( corresponding to 
photoemission ) the corresponding hermitian conjugate 
operators should be taken. In real space using Eq.(2) 
one obtains: 



di = D x c ia {\ - n i9 ) + D 2 c i(y n la 

a 1 = P\Cia(l - Ilia) + P 2 Cian is 



Transforming Fourier and using Eq.(3), the operators en- 
tering Eq.(6) are defined. As will be discussed in section 
IV, the asymmetry of Eq.(7) and also the Hamiltonian 
Eq.(5) under electron-hole transformation leads to a dif- 
ferent behavior of the spectral densities under electron or 
hole doping. This difference is apparent in the cuprates 
and in the three-band Hubbard model results, but it is 
absent in the ordinary one-band Hubbard model jp 
The optical conductivity can be written as 1119,111 




<j(lS) = 2ttD5(uj) + o reg (uj) 



(9) 



where 



<W") = ^- \(n\j x \g)\ 2 S(E n -E g -w) (10) 



V is the volume of the system, j x the x component of the 
current operator. The Drude weight D can be calculated 
in two ways: either using the mean value of the kinetic 
energy in a given direction (T x ): 



D 



-T 



2V 



v 2-? E n - 



E„ 



(11) 



or from the second derivative of the energy with respect 
to the vector potential A in a system with periodic ( 
or arbitrary as described in the next section ) boundary 
conditions 



D = 



d 2 E(A) 



2V dA 2 



A=0 



(12) 



The current operator can be derived in a standard way, 
from dH/dA, when each hopping term c^ ia Cj a is multi- 



(7) plied by a phase factor F,xp[ieA(xi —Xj)/{cK)] |l6|Jl7[]. In 
the case of the effective Hamiltonian Eq.(5), the current 
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can be divided in three terms proportional to the three infinite square lattice, the allowed total wave vectors K 
different types of hopping terms: should satisfy 



Jx = Jaa + Jab + Jbb 



(13) 



with 



AT0i + 2ni7r = K Vi 
N(/)2 + 2n 2 7r = K ■ V 2 



(15) 



3Aa{t) = ~^ f AA y^(4 g Cr+a;g ~ 4+xct c W)( 1 - n r ^)(l - n r+x . s ) (14) 

<7 

jBB(r) = —frtBB 2_J,c r <j c r+xo - c;, +xa c rcr )n r ^n r+XtS 



An alternative derivation of the current operator for 
the effective one-band model is to obtain first the current 
operator of the multiband model p7j , then project the 
result into the low-energy subspace and map it onto the 
one-band representation. Both derivations lead to identi- 
cal results if t 



pp 



0. For realistic values of t pp , there are 



small differences between the resulting coefficients enter- 
ing the current jaa, Jab and Jbb- For example while 
for the set 2 of parameters of (see table I ) we obtain 
tAA = tAB = —0.38, tgg = —0.33, in the second deriva- 
tion these values should be replaced by -0.32, -0.35 and 
-0.36 in the respective currents. This difference is prob- 
ably due to the different way in which excited ( mainly 
local triplet ) states are taken into account in lowest or- 
der in both procedures when t pp ^ 0,^6| . In this work 
we use Eq.(14). 



III. NUMERICAL METHODS AND AVERAGE 
OVER BOUNDARY CONDITIONS 

We have evaluated Eqs.(6), (10) and (11) using the 
now standard continued-fraction expansion of these equa- 
tions with the Lanczos method J62| , in a square cluster 
of 10 unit cells, integrating the result over the different 
boundary conditions jll],|4| . Particular boundary condi- 
tions are specified by two phases (0i,0 2 ), < 0j < 2n 
in the following way: an infinite square lattice is di- 
vided into square 10-site clusters, and the sites which 
are at distances nVi + mV 2 with n, m integers and 
Vi = (3,1), V 2 = (—1,3) are considered equiva- 
lent ( see Fig.l). After choosing a particular 10-site 
cluster, each "time" the hopping term makes a partic- 
ular jump out of the cluster, it is mapped back into 
the cluster through a translation in one of the vectors 
— Vi , Vi , — V 2 , or V2 and the wave function is mul- 



tiplied by e^ 1 , e 



-i<t>2 



respectively. When 



(0i,02) = (0,0) this is equivalent to the usual periodic 
boundary conditions. It is easy to see that when we im- 
pose the wave function for N particles to be an irre- 
ducible representation of the group of translations of the 



where m, n 2 are integers. Solving for K x , K y one obtains: 



10 v v v ' 10 



(16) 



When (0i,02) = (0,0), varying the integers ni 7 n 2 , 
Eqs.(16) give the 10 inequivalent wave vectors allowed by 
periodic boundary conditions. When (0i,02) is allowed 
to vary continuously the whole reciprocal space can be 
swept. This allows us to obtain more information from 
the finite-size cluster calculations. Also, while the spec- 
tral densities and optical conductivities being calculated 
are given by a sum of delta functions for each (0i, 02) ( 
Eqs.(6) and (9)), integrating the result over (0i, 2 ) leads 
to continuous functions. In practice we have replaced the 
integral by an average over a square mesh of up to 8x8 
points ( up to 34 non-equivalent by symmetry for the op- 
tical conductivity) in 0i , 2 space until a fairly smooth 
function was obtained. 

However, the procedure described above has some 
shortcomings when applied to the optical conductivity. 
For general boundary conditions there is a spontaneous 
current in the ground state, invalidating the derivation 
of Eqs.(9) to (12) |l^,|l^]. (Such a current can not exist 
in the thermodynamic limit, even if it were allowed by 
symmetry |6^j6^1 ). In finite systems, this current is zero 
for periodic boundary conditions at "closed-shell" fillings, 
(including half- filling in our cluster), and for half- filling 
and any boundary conditions in the U — > 00 limit (as 
in the t — J model |24[] ), since in this limit the charge 
dynamics is suppressed. Our point of view is that the 
average of Eqs.(10) and (11) over boundary conditions 
in a finite system, is an approximation to the correct 
expression in the thermodynamic limit, which converges 
faster with system size than the result of Eqs.(10) and 
(11) or (12) for fixed boundary conditions. Since we 
have not done finite-size scaling in two dimensions, we 
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can not prove this statement. However, it is supported 
by our one-dimensional check summarized below. Fur- 
thermore, as we shall see in section V, the results us- 
ing the averaging procedure look reasonable, while the 
Drude weight obtained in larger systems using periodic 
boundary conditions show unphysical negative values at 
half- filling |2^,^]. Previous calculations of the optical 
conductivity of the t — J model using twisted boundary 
conditions also support our procedure p4j . 

We have checked the method of averaging over bound- 
ary conditions applying it to the optical conductivity 
a(<jj) of one-dimensional Hubbard rings, for which finite- 
size scaling can be done, and some exact results, like the 
charge gap, are available p3| . The results for U/t = 4, 6 
and 8 and several system sizes L are represented in 
Figs. 2-4 respectively. A funny feature of these curves 
is the presence of oscillations in a(uj). A comparison of 
the result for rings of different sizes shows that this is a 
finite-size effect. Calculations using spin- wave theory in 
the strong coupling limit give a smooth <j{lo) |p6| . For- 
tunately no signs of these oscillations are present in the 
two-dimensional case. By comparison, the strong cou- 
pling limit of the Hubbard model, and the t — J model 
with one hole in systems as large as 19 sites with periodic 
boundary conditions shows a few delta functions p6[ . 

The presence of a peak at zero frequency is a con- 
sequence of the already discussed spontaneous current 
in the ground state for general boundary conditions. 
Clearly, this is also a finite-size effect and the height of 
the peak decreases with system size and also with in- 
creasing U/t. Thus, we have suppressed this peak in the 
two-dimensional results of section V. 

We have also noticed that the gap in <j(oj) converges 
much faster to the thermodynamic limit E g than the cor- 
responding result using fixed boundary conditions. In ta- 
ble III, we compare the exact gap, given by the integral 

©; 

%L = ™ f dy (17) 

t U A sinh(^) ^ [ 

against its value ( taken as the position of the lowest 
frequency peak in the spectra of Figs. 2-4 ) for a ring 
of 8 sites after averaging over boundary conditions and 
at fixed antiperiodic boundary condition. Clearly, the 
averaging procedure speeds up the convergence toward 
the Loo limit. 

IV. THE SPECTRAL DENSITIES 

In Fig. 5, we show the Cu and O spectral densities p(u) 
in the hole representation ( photoemission and inverse 
photoemission correspond to u> > €f and to < 6f respec- 
tively, where is the Fermi energy ) for several doping 
levels and the set of parameters 1 ( see Tables I and II 



). These spectral densities have been calculated as de- 
scribed in section III, using Eq.(6) with Ok replaced by 
d\, pfr for Cu and O photoemission, and by dk, P% for 
inverse photoemission. The result for O has been multi- 
plied by a factor 2 to represent the total contribution of 
both 2p a orbitals of the unit cell. Also shown in Fig. 5 are 
the extended Hubbard spectral densities corresponding 
to the operators c\.(u> > ep) and Ck(oJ < £f)- Note that 
in constrast to the usual Hubbard model, the extended 
Hubbard spectral densities are not symmetric under hole 
(x > 0) or electron (x < 0) doping. This is due mainly 
to the fact that tAA ^ tsB- However, the asymmetry 
of the Cu and O spectral densities is much stronger. In 
agreement with experiment, for hole doping, the inverse 
photoemission spectrum near ep shows a larger contri- 
bution of O states than Cu d states, while for electron 
doping, the photoemission spectrum near the Fermi level 
is dominated by Cu states. 

A pseudogap in the spectral densities persists with 
doping and allows to separate the electronic states in two 
bands at lower and higher energies than this pseudogap. 
The lower band corresponds to states with low double 
occupancy and therefore, the spectral densities are dom- 
inated by the contribution proportional to Ci CT (l — nig-) 
in Eq.(7); thus the Cu and O spectral densities in the 
lower band first approximation proportional to 

the respective coefficients( B\ and B 2l see Eq.(8) ) in the 
one-particle ground state of the cell. However, the states 
of the lower band contain some admixture of double oc- 
cupied states, and this introduces an energy dependence 
of the relative weight of Cu and O states. The states of 
lower energy (to ~ —1) have a larger Cu content than 
those of the rest of the lower band. Similarly, as dis- 
cussed previously by Feiner [Q , the Cu and O contents 
of the upper band are mainly determined by D 2 and Pi 
respectively ( see Eq.(8)), which depend on the structure 
of the Zhang-Rice singlet \i2) and the ground state of the 
cell with one hole \ia) (see Eq.(4) ). 

From Eq.(6) applied to the Hubbard operators Ck and 
c\, it is easy to obtain the total inverse photoemission 
N~ and photoemission N + spectral weight: 

p c (u)du) =< m\ > + < Tin >= 1 + x ( 18 ) 

-00 

/>oo 

N + — / p c +(ui)doj — 2— < rii^ > — < nn >= 1 — x 

These sum rules are the same as those corresponding to 
the usual Hubbard model. However, using Eqs.(6) and 
(7) the following sum rules for the Cu orbital d x 2_ y 2 and 
" antibonding" O Wannier function a are obtained: 

N- = (n d i + n di ) = (1 + x)D\ + 2d(D 2 2 - D\) (19) 

N+ = 2{d - x)D\ + (1 + x - 2d)D\ 

N~ = (n aT + n al ) = (1 + x)Pl + 2d(Pl - Pf) 

N+ = 2{d - x)Pl + (1 + x - 2d)Pl 
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where d =< n^riii Adding the four functions, the 
total spectral weight becomes: 

N t = {l+x)(l + Dl + P$)-2x (20) 

Since D^ + P^ < 1, then N t < 2. The rest of the spectral 
weight has been projected out of the low-energy effective 
extended Hubbard model and it is contained at higher 
energies. The spectral weights given by Eqs.(19) as a 
function of doping are represented in Fig. 6 for a typical 
set of parameters of the three-band model. In Fig. 7, we 
show the coefficients Pi , Di of the effective operators ( 
Eq.(8) ) as a function of the three-band model parameter 
A. We also show the function [I — (D\ + -P 2 2 )L related by 
Eq.(20) to the amount of spectral weight projected out 
of the effective low-energy Hamiltonian. From the slopes 
of and N~ as a function of doping ( Fig. 6 ) one 
realizes that when doping the stoichiometric compound 
with electrons, the latter occupy mainly Cu states, while 
for hole doping, holes occupy mainly O states, and for the 
parameters of Fig. 6, part of the Cu holes are transferred 
to O holes. This is a well known effect of the Cu-0 
repulsion U p d |3jJ. 

From the above discussion, it is clear that the main fea- 
tures ( peaks and valleys) of the Cu and O spectral den- 
sities are present in the corresponding extended Hubbard 
result ( obtained with the Ck and c k operators ) , and as a 
first approximation, the former can be obtained from the 
latter modifying the spectral weights of the lower and up- 
per Hubbard bands with the coefficients given by Eqs.(8) 
and represented in Fig. 7. Thus in the following we discuss 
the general properties of the extended Hubbard spectral 
density ( independently of the mapping procedure ) and 
the spectral weight of both extended Hubbard bands. 

In Fig. 8 we show the evolution with hole doping of the 
spectral density of states of the extended Hubbard model. 
As discussed previously in the case of the ordinary Hub- 
bard model |P2[ , after doping, states appear in the gap of 
the stoichiometric compound and the Fermi level jumps 
to the upper Hubbard band. According to our results, 
these states are distributed evenly inside the gap and 
there is not an impurity band at a defined energy. Also, 
the pseudogap persits for all dopings, and both bands 
evolve smoothly with doping, with a noticeable transfer 
of spectral weight from the lower to the upper band. In 
Fig. 9, we show the same results for the set of parameters 
3 of Tables I and II, for which the ratio of the on-site 
Coulomb repulsion to hopping parameters is larger. In 
this case, the weight of the states which appear in the 
gap after doping is much smaller, and also the amount 
of spectral weight transfer between the bands is smaller. 
Note also and by contrast to the usual Hubbard model 
the strong asymmetry of p(uj) at x — due to tAA ^^bb- 

The transfer of spectral weight in the Hubbard model 
has been discussed previously |j],[30) and is important in 
the analysis of x-ray absorption spectra in the cuprates 



In the strong-coupling limit ( large U ), or in the 
Hubbard III approximation, a static spectral weight of 2x 
is easily obtained ( 1 — x states lie in the upper Hubbard 
band). However, as the hopping increases, a positive 
dynamical contribution to the spectral weight becomes 
important. In Fig. 10 we show this dynamical contribu- 
tion for the extended Hubbard model and the three sets 
of parameters of Table II. It is apparent that for realistic 
parameters for the cuprates, the dynamical contribution 
is very important, and as could be seen comparing Figs. 
8 and 9, it is more important for smaller ratios U/t m , 
where t m = {tAA + ^ab + tss)/4 is the average effective 
hopping. 

V. THE OPTICAL CONDUCTIVITY 

For the evaluation of Eqs.(10) and (11), we have taken 
the distance between planes d± = 6.64A, corresponding 
to La2Cu04. Thus, the Drude weight D and the regular 
part of the optical conductivity <j reg (uj) are proportional 
to the constant 2e 2 /{M 1 _) = 3.7 x lO^^cm" 1 . The 
resulting Drude weight and average value of the kinetic 
energy as a function of hole doping are represented in 
Fig. 11 for the choice of parameters 3 of Table II. The cor- 
responding results for the optical conductivity are shown 
in Fig. 12. The average over all boundary conditions has 
been replaced by a discrete sum over a square mesh of 
8x8 points (0i,</>2) ( 34 non-equivalent by symmetry ). 
For more than 16 points the results are practically inde- 
pendent on the number of different boundary conditions 
taken. 

For the semiconducting system (x = 0), there is a spu- 
rious peak at u> ~ 0.4. This is a finite-size effect related 
with the fact discussed in section III, that the ground 
state for general boundary conditions carry a current for 
finite U . The peak corresponds to transitions to excited 
states which differ from the ground state mainly in the 
spin arrangement. For realistic and large U the magni- 
tude of the spurious peak is small. In Fig. 13, we show the 
optical conductivity for a smaller value of the effective U 
after averaging over 32 (</>i, fa) points ( 18 non-equivalent 
by symmetry) . The magnitude of the spurious peak in- 
creases, but all the other features are very similar be- 
tween them, and qualitatively similar to previous results 
obtained in periodic 4x4 clusters Ej , However, our pro- 
cedure of averaging (Eq.(ll) ) over boundary conditions 
leads to more reasonable ( small and positive ) values 
of the Drude weight for the insulating system: 5 x 10~ 5 
and 3 x 10 _3 f2 _1 cm _1 for the set 3 and 2 of parameters 
respectively. 

For x < 0.2, the results are in semiqualitative agree- 
ment with experiments J^Q- Excluding the spurious 
peak near 0.4, the optical conductivity of the insulat- 
ing compound shows a gap of ~ 2eV and a character- 
istic shape at larger energies. Moreover, if the scale 
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of energies is modified by a ~ 3/4, the main features 
of Fig. 13 for x — coincide with the observed spec- 
tral dependence of the persistent photoconductivity in 
YBa 2 Cu306.38 |@- As the system is doped the spec- 
tral intensity above ~ 2eV decreases and at the same 
time, a mid-infrared peak appears. Although this low- 
frequency peak has been discussed before, its nature has 
not been clarified yet. Recently, an explanation based on 
the string picture of the t — J model has been proposed 
p9[ . In order to shed light on the origin of the different 
contributions to the optical conductivity, we have sepa- 
rated the contributions of the three currents, jaa, 3bb 
and jAB, and compared them with the spectral density, 
as shown in Fig. 14. For hole-doped systems, the contri- 
bution of jAA is negligible. From the analysis of Fig. 14, 
we conclude that the optical conductivity at energies of 
the order of the gap and above corresponds to transitions 
between the lower and upper extended Hubbard bands, 
originated by tAB ( which corresponds to a hopping be- 
tween two nearest-neighbor singly-occupied sites ). In- 
stead, the mid-infrared peak corresponds to one-particle 
excitations inside the upper Hubbard band ( from below 
to above the Fermi level ) , originated by ts b or in other 
words ( in agreement with Ref. p9| ) to movements of 
the added holes without changing the amount of doubly 
occupied sites. The difference between the total optical 
conductivity and the contributions of jAB and jbb rep- 
resented in Fig. 14 is mainly due to cross terms involving 
both currents. 

In Fig. 15 we show the frequency dependent effective 
number of carriers defined by 

N ef (uj) = ^^ (" o{J)du', (21) 

for the two sets of parameter used before in this sec- 
tion. In order to compare with experiment, we took the 
lattice parameter of the CuC>2 planes a — 3.78A, which 
corresponds to La2CuC>4. As well as the results shown 
in Figs. 12 and 13, the agreement with experiment is 
good for x < 0.2. Experimentally, for x ~ 0.3 there 
seems to be an abrupt change of regime to an uncorre- 
cted metal, with the Drude peak as the only significant 
feature of the optical conductivity. Instead, our results 
do not show such a transition, but only a smooth evolu- 
tion. Possible reasons of this discrepancy are discussed 
in the next section. 

VI. SUMMARY AND DISCUSSION 

Using an effective extended Hubbard model for the 
cuprates, we have calculated the Cu and O spectral den- 
sity and optical conductivity in a square cluster of 10 
unit cells, averaging over all possible boundary condi- 
tions. To our knowledge, this is the larger cluster for 



which Cu and O spectral densities are calculated using 
the Lanczos method. The averaging procedure allows to 
obtain continuous spectral densities and optical conduc- 
tivity, and leads to more reasonable values of the Drude 
weight near stoichiomctry. 

We obtain that the states which appear in the gap after 
doping are distributed evenly inside it, without building 
a new band of midgap or impurity like states. We also 
obtain that a marked pseudogap persists for all dopings, 
separating two bands. When the insulator is doped, the 
Fermi level jumps into the band which corresponds ac- 
cording to the sign of the doping. As doping proceeds, 
spectral weight is transferred to the band which contains 
the Fermi level from the other one. An analysis of the 
different components of the current operator allows us to 
conclude that the mid-infrared peak in the optical con- 
ductivity observed in doped systems is related to intra- 
band transitions across the Fermi level. The shape of 
the Cu and O spectral densities is consistent with the 
fact that for hole doping, holes enter mainly O p a states, 
while Cu states are occupied on electron doping. 

In general, these results are in agreement with experi- 
ment. There is however a disagreement with some optical 
experiments [||,[)| which indicated almost no shift of the 
Fermi level with doping, as mentioned in section I. Wc 
should mention that, as discussed by Hybertsen et al. 
the substitution of La +3 by Sr +2 creates a potential 
which might localize the doped holes by small doping, 
creating an impurity band of acceptor levels like in or- 
dinary semiconductors. To include these effects it would 
be necessary to generalize the three-band Hubbard model 
and the effective models derived from it. 

We do not obtain a deep reconstruction of the elec- 
tronic structure around hole doping for x ~ 0.3 as con- 
ductivity j|, Hall and other experiments suggest. A 
closing of the Mott-Hubbard gap was expected for sim- 
ilar dopings [0. This closing in fact can occur in the 
three-band Hubbard model and the effective extended 
one band model for U p d ~ 3 — 4. For these large values 
of Upd the effective U of the one band model decreases 
with doping and a large transfer of holes from Cu to O 
takes place J35|. However, these values of U p d are too 
large compared to the most accepted ones. For realistic 
parameters of the three-band model, the effective model 
parameters are independent of doping |35|| . This does 
not exclude the possibility that the three-band parame- 
ters change near x ~ 0.3, as a consequence of screening 
caused by the added holes, for example. Nevertheless, 
the fact that there are always two bands as a function 
of doping is in agreement with x-ray absorption experi- 
ments in La2- £C Sr :! ;Cu04 [BUT]]. 

Due to the small size of the cluster, we have looked 
neither for the presence of a Kondo-like peak near the 
Fermi energy p5[ , nor the effects of excitons in the opti- 
cal conductivity, which are present in the effective model 
for realistic and large Vij pi. A signal of the presence 
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of these excitons is a smaller gap in the optical conduc- 
tivity than in the spectral density. For example, for zero 
hoppings, it is easy to see that the gap in the spectral 
density is given by 

E gap = U + 8(V 12 -Vu) (22) 

However, if the hole and doubly occupied sites are near 
each other, the excitation energy is 2Vi2 — V\\ smaller. 

We conclude that for the most interesting doping lev- 
els and the most accepted values of the parameters of 
the three-band model, the general features of the spec- 
tral properties of this model and of the cuprates, can 
be described well using an effective extended one-band 
Hubbard model. 
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FIGURE 1: 

Scheme of the supercell ( cluster ) containing 10 unit 
cells used in our calculations. The vectors Vi and V 2 
generate the whole lattice by succesive translations of 
the supercell. 

FIGURE 2: 

Average optical conductivity as a function of frequency 
for Hubbard rings of different length and U/t = 4. 

FIGURE 3: 

Same as Fig.2, for U/t = 6. 
FIGURE 4: 

Same as Fig.2, for U/t = 8. 
FIGURE 5: 

Inverse photoemission (w < ep) and photoemission (w > 
ep) hole spectral densities for Cu ( dashed line), O ( 
full line ) and Hubbard operator ( dotted-dash line ) , for 
different hole dopings x. The vertical line indicates the 
position of the Fermi level ep. Parameters are given by 
the set 1 of table I. The average over boundary conditions 
have been replaced by a discrete sum over 4x4 (cfrijfo) 
points ( 6 independent points). 

FIGURE 6: 

Total photoemission (N + ) and inverse photoemission 
(N~) spectral weight for Cu ( ) and O (N p ) as a 
function of hole doping. Parameters are Ud = 9, A = 
■i,U p d = l,i PP = 0.5,U P = 4. 

FIGURE 7: 

Parameters which define the effective low-energy Cu and 
O operators according to Eq.(8) as a function of A. Also 
shown is 1 — D\ — P|, which is related with the high- 
energy spectral weight projected out of the low-energy 
effective Hamiltonian ( sec Eq.(20)). Other parameters 
as in Fig. 6. 

FIGURE 8: 

Extended Hubbard inverse photoemission (u> < ep) and 
photoemission (oj > ep) hole spectral densities for dif- 
ferent hole dopings x. The vertical line indicates the 
position of the Fermi level ep. Parameters are given by 
the set 1 of table II. 

FIGURE 9: 

Same as Fig. 8 for the set 3 of table II. 
FIGURE 10: 

Total amount of occupied hole states in the upper band 
W minus two times the hole doping function of 

x for the three sets of parameters of Table I and II. For 
U -» oo, W - 2x = 0. 



FIGURE 11: 

Drude weight and average kinetic energy of the effective 
model Eq.(5) as a function of doping. Parameters are 
given by the set 3 of table II. 

FIGURE 12: 

Optical conductivity for different doping levels and the 
set 3 of parameters of table II. The delta function con- 
tribution at u = has been replaced by a Lorentzian of 
with 0.01 to facilitate the comparison with experiment. 

FIGURE 13: 

Same as Fig. 12 for the set 2. 
FIGURE 14: 

Bottom: spectral densities of I mag « >> ( con- 
tributions to the optical conductivity ) of the partial cur- 
rents Jab, Jbb and the total result ( see Eqs.(10) and 
(13) ) for x = 0.2 and the set 2 of parameters of Table 
II. Top: corresponding hole spectral density. The arrows 
indicate the transitions which give rise to the different 
contributions to the optical conductivity. 

FIGURE 15: 

Effective number of carriers defined by Eq.(21) as a func- 
tion of frequency for different doping levels and the sets 
of parameters 3 and 2 of Table II. 
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TABLES 

TABLE I : Different sets of parameters of the three-band 
Hubbard model H^b ( Eq.l) in units of t p d- 



Set 


A 


u d 


Upd 


tpp 


u p 


1 


2.0 


7 


1 


0.5 


4 


2 


2.5 


9 


1 


0.5 


4 


3 


4.0 


7 


1 


0.5 


4 



TABLE II: Effective one-band parameters for the three 
sets of parameters of H 3 b listed in Table I. The last col- 
umn gives the ratio of the effective on-site repulsion U to 
the average effective hopping t m — (tAA + 2tAB +tss)/4. 



Set 


U 


tAA 


tAB 


tBB 


Vii 


Vl2 


V 22 


U/tm 


1 


2.15 


— 0.382 


-0.385 


— 0.337 


0.123 


0.170 


0.178 


5.78 


2 


3.01 


— 0.377 


-0.324 


— 0.338 


0.119 


0.162 


0.181 


8.85 


3 


3.37 


-0.268 


-0.327 


-0.337 


0.066 


0.125 


0.156 


10.71 



TABLE III: Exact gap E g compared with the gap E a g v 
obtained after averaging over boundary conditions and 
for fixed ( antiperiodic ) boundary condition for a 
ring of 8 sites and different values oiU/t. 



u/t 


E 9 


£-(8) 


^(8) 


4 


1.2868 


1.98 


3.68 


6 


2.8928 


3.50 


4.90 


8 


4.6796 


5.21 


6.47 
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